Semi-leptonic B-decays and the two-pion distribution amplitudes 
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Abstract. We show that the semi-leptonic decay B + — > -k + -k~ l + i>i can be used as a source of information 
for two-pion distribution amplitudes. The connection between these amplitudes and the B-meson decay 
width is achieved by the light cone sum rule method. We show the relevant distribution amplitudes and 
give the formula for the decay width. 

PACS. 12.38.Lg Other non-perturbative calculations - 13.60.Le Meson production 



1 Introduction 

Recently, two-pion distribution amplitudes (2-7rDAs) have 
received considerable interest [1,2] because of their rela- 
tion to skewed parton distributions [3]. From an exper- 
imental point of view the 2-7rDAs have been mostly dis- 
cussed for the reaction 77* — > 7T7t [4,5,6,7], where the 
factorization has been proven in Refs. [8,9], and also in 
hard exclusive electroproduction [10,11]. Here, we want 
to add another type of reaction which could provide valu- 
able insight in the dynamics of the formation of two pions, 
namely semi-leptonic B-decays. 

The reaction B — > irwlv is an alternative source of in- 
formation about the 2-7rDAs as compared to 77* — * ttit 
because new structures arise due to the fact that the semi- 
leptonic weak decay induces an axial vector current in ad- 
dition to the vector current and the fact that the B-meson 
is a pseudoscalar particle. 

The method used here to connect the 2-7rDAs with the 
B-meson decay width is the method of Light Cone Sum 
Rules (LCSR) where the decay amplitude B — ► rnviu is 
related to the light-cone OPE of the corresponding cor- 
relation function where factorization into the 27rDAs and 
the hard amplitude is guaranteed by the large virtuality 
of the off-shell currents. The light-cone sum rule method 
applied here is essentially the same as used in B — > -niu 
[12,13] and B — > piv [14], with the only distinction that 
the 27rDAs enter. The advance of B-factories may yield 
a lot of new experimental data on the decay B — > ttit£i>, 
where explicit models for the various 2-7rDAs entering in 
this process may be tested and could yield a deepening 
understanding of the non-perturbative multi-particle dy- 
namics which lies behind these generalized distribution 
amplitudes. 



2 The method 

2.1 Kinematics 

We consider the process B + — > 7r + 7r~i? + tv- The kinematics 
of the process is given by B(q) — > £(p e ) + v{q' — Pe) + 
7r + (fci) + 7r~(fc2) and there exist two light-like vectors n + 
and n~ with n + n~ — 1/2 such that: 

q = (m B , 0,0,0) 

/ 1— 

q r 



rriB [n 



9 



/+ - 
q n 



(1) 



in the rest frame of the B-meson. We can now make a 
division into 'good' (+) and 'bad' (-) components, where 
the 'bad' components can be neglected requiring ms S> 
q'_ > and m| > P 2 = (fci + k 2 ) 2 = W 2 . This require- 
ment is necessary to ensure the factorization in the ap- 
proach of the LCSR technique at the stage where the vir- 
tual amplitude is factorized in a hard scattering part and 
the 27rDAs. The factorization follows in complete analogy 
from the 77* case. Under these circumstances P 2 /m' B 2 , 
using m' B = ttib — q' _ can be considered as a small ex- 
pansion parameter and we can decompose: 



p2 

P = q - q' = ki + k 2 = m' B n + H r n 

m B 

-P 2 

ki = Qm' R n + + ( — — n + K± 
m B 



,; 2 = (m' B n + + C^-n - K± , 



P 2 



(2) 



using £ = 1 — £. In this way we have set up a similar light 
cone decomposition as in the case 77* — > tt-k. 

For the light cone sum rule technique to apply we need 
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Fig. 1. Diagrammatic representation of the correlation func- 
tion. The wavy line represents the external 'B-meson current' 
Jb and the dashed line the weak current J wcak . 



to consider a situation where the B meson is off shell. 
One can achieve this in the frame- work of the kinematics 
discussed so far by simply changing: 



niBn^ H n 

m B 



(3) 



The two light-like vectors n + and n are still the same as 
in Eq. (1). 



2.2 The correlator and the distribution amplitudes 

For the application of the LCSR approach one considers a 
correlation function of a pseudoscalar and a weak current 
with a quark content that corresponds to the one of a B ± 
meson: 



g 

twist — 3 : 



Jo 



#W = W [ dze™(zX,W 2 )e 
Jo 

^ a u J 



2\AzP-x 



(6) 



using = fc 2M — fci M . Especially, we find P ■ R = and 
R 2 = 4m 2 . — W 2 . We only need one of the two Lorentz 
structures 10^75 and a^ v , as they are related to each 
other by: 

(7) 



a = o £ CTa/375 • 



With this decomposition we obtain for the correlation 
function Eq. (4): 



Tfj, — ie^pRq' V u b 



+ 



+ 



2m b g^(z,C 1 W 2 ) 
((Pz + q' f -m 2 Y 

f:i(z,(,W 2 )/(4Wf 
(Pz + q') 2 - m 2 

Vub 



m b PJ™(z,(,W') 



(P z + q') - m 2 

(Pz + q') u We^{zX,W 2 ) 

P»R ■ q' - R^P ■ (Pz + q') 
AW 



JZl(zX,w 2 ) 



.(8) 



T fl (q,q')=i J d 4 xe-^ x (0\T[J B (0),Jj: eak (x)}\n + n-(P)) 
Jb{x) = u(x)ij 5 b(x) 



JJ^(x)=b(x)i^(l- l5 )u(x) 



(4) 



In this correlation function the B-meson is interpolated 
by a current Jb with four-momentum q, which is off shell, 
jweak j g wea k jj _^ u transition current. The diagram- 
matic representation of T M is shown in Fig. 1. Depending 
on their angular momentum the two pions in the final 
state can have odd or even parity. With the definition: 



(n + TT-(P)\q(x)rq(0)\0) x 2 =0 =$\p , 



(5) 



we can parameterize the structures that will occur in our 
correlation function to leading-twist accuracy: 



2.3 Calculation of the decay amplitude in the LCSR 
approach 

The next step is to relate the decay amplitude to the distri- 
bution amplitudes described to far. For this purpose one 
inserts a complete set of states with B-meson quantum 
numbers between the currents [12,13]: 



F^q 2 ) = J d 4 xe 



-%q x 



(0|a(0)z 75 fr(0)|g) 
q 2 - m% 



x(B\b(x)i~f^(\ - 7 5 )w(x)|7T7r) + . . . 

-M u + ... 



fB.m 2 B 



m b (q 2 - m%y 



(9) 



using: 



twist — 2 : 



g 



2W jo 



dzfr q (z,CW z )e 



2\izP-x 



f 

Jo 



dz fr q (z,CW 2 )e 



2\izP-x 



e„ kl k 2 x I dz g^(z,(,W 2 )e 
'0 



2\izP-x 



M u = / d A xe- l i x (B\b(x)j^(l - l5 )u(x)\irir) . (10) 



The ellipsis indicates all the other hadronic states which 
in the end we will suppress by Borel transformation. The 
next step is to define the discontinuity Disc: 



Disc [T^(q 2 



1 

2ni 



(T^q 2 -ie)-T^q 2 +ie)) , (11) 



Martin Maul: Semi-leptonic B-decays and the two-pion distribution amplitudes 



3 



and isolate using the standard duality approximation: 

pSQ s—m?" 

\ d S Disc[T,( S )]e-^-^-=M,. (12) 

J ml ™ B JB 

Now one can expand M.^ into an orthogonal-system: 
M-n = J rf 4 x e - 4 "' 3; (B|6(x)7 Al (l- 75 ) U (x)|7r7r) . 



■Mfi = M X P„ + M 2 R„ + M 3 <# 

n j R-q' P-q' 
it,, 



2W 



% = % 



R 2 



P, ■ 



pi v- 



(13) 



With this decomposition one gets the following LCSR re- 
sults: 



Mi 



-L 



1 dz 



M> 



-L 



z m%f B 

-lw( s + {2z - l)w2 - q ' 2 ) e ™ 

i , - a ~< 
dz e m 2 nib 

2 



9[c(z,af)] 



z m B .fB 



R 2 



-±^^+(2z-l)W 2 - q ' 2 )r u l(z) 



0[c(z,s*)} 



M, 



Ma = 



1 dz 



z m B h 



WeZ*(z)e[c(z, S *)} 



1 dz e nib 
o z ™? B f B 



™fZl{z)9[c{z,s*)] 



4m b W 2 9 r (z) I ^0[c(z, S f)] + 6[c(z, s*)} 



-6[-c{z, ml)] 



(14) 



using: 



2(R.q')(s) = ((-0[m B q' + - 



sq' 



m B 



zzW 2 +ml — zq^ 

z l 

c(z, sf ) = zs B - m 2 b + zq' 2 - zzW 2 . (15) 

and 8 functions arise from the continuum subtraction 
[14]. More precisely we use: 



-i T r , r s ° ds'e — ^f(z) 

— lm dz — - -r: '-^ 

Jo J m l (Pz + q') 2 -m 2 + ie 



1 dz 



eST- f(z)G[c(z, s*)]G[-c(z, m 2 )] 



-1 f 1 [■ 
Im I dz I 

K Jo Jrr, 



d_ 

da 

L 



Im I dz i 

n n Jo Jn 

a=0 



'<> ds'e ^~f(z) 
2 [(Pz + q') 2 -m 2 +ie} 2 

s' -m 2 

ds'e ^ S 'f{z) 



1 dz 



= — i — e m j 
)o z 



'f(z) 



zM 2 



2 (Pz + q') 2 - m 2 b + ie + a 

e[c(z,8f)]+S[c(z, 8 f)] 



-5[-c(z,ml)] 



(16) 



Note that 0[-c(z,m 2 )] = 6>[1 - z], and 5[-c(z, m 2 )] = 
6(1 — z)/(m 2 + zW 2 — q' 2 ). In order to obtain the square 
of the decay amplitude, we have to multiply with 
the leptonic scattering tensor: 

Lftft' = ^u s (p e )7 M (l ~^)u s >(p v ) 

ss' 

x u s '(p i ,)7 (U '(l - -f 5 )u s (Pe) 

= 8 [{q'f, - Pe»)Pen> + (q'f,> - Pep,')Pep, 

- 9w'(j>e ■ q' - m 2 e ) +ie wVpe ] . (17) 
Then one gets for the matrix element: 



IM 



2 _ m^m;, 
~ W - q' 2 ) 2 



16(MiP • Pe + M 2 R ■ Pe + M 3 q" ■ p e ) 

X (MiP • q' + M 2 R ■ q' + M 3 q" ■ q') 
-16(MiP • Pe + M 2 R ■ Pe + M 3 q" ■ p e f 

+4M 4 2 ^ 
4 w 2 

xa" 2 (m 2 {R - Pe)2 (P - P ^ ( q "- p ^ 



R 2 



W 2 



-8 ^M{P 2 + M 2 R z +M 2 q" 2 - Mr 
x(p e ■ q' - m 2 e ) 

-8^e^, qlp ^ PRq/ (MiP^ + M 2 R^ + M 3 q'^) 
1 

X (M 2 W - q*y ■ 



-,112 



f 2 R „//2 



AW 2 



(18) 



Here we have dropped the weak coupling (g/(2v / 2)) 4 for 
simplicity. We will add it later to the phase-space element 
in Eq. (23). The full expression Eq. (18) is rather com- 
plicated. To simplify it we make first use of the fact that 
P 2 /m' B 2 is small, so that we can throw away all 'bad' 
components connected with n~ which corresponds to an 
expansion in P 2 /m' B 2 . Then we get approximately the fol- 
lowing expressions: 
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2P 


■q' 


P 2 

= m' B q' H r q'~ - 


-> m B q 






■q' 


= {q' 2 + ml) 






2P 


Pe 


p2 

= m BPi + —TPe - 


m B p+ 




2R 


■q' 


= (C-C) (rn' B q' + - 


P 2 ,-\ 

—q 

m B J 


-> (1 - 2C)m B(J ' + 


2R 


Pe 


= (C-C) (rn' B Pe + - 


p2 X 
—Pe 

m' B J 


+ 4K ± • p e± 



(1 - 2()m BPe + + 4K ± • p e± . 



(19) 



Here we have defined pf = 2p e ■ n + . The independent 
variables are now q' + , p+ and p e ±. and one should note 
that 

K ± 2 = aW 2 -ml . 
We can define the dimensionless quantities: 



(20) 



m B 

Asmfl> q'~ we find also: 



y 



Pi 



<C x . 



(21) 



(22) 



To simplify the expression Eq. (18) further, we consider 
the limit where q' becomes quasi light-like, i.e. q' 2 = 0. 
Then the leptonic part of the decay process simply factor- 
izes in analogy to the Weizsacker Williams approximation 
in photoproduction, see e.g. Ref. [15]. In this limit q' is 
quasi collincar to p e , and one obtains after integration over 
the angle related to p e ± and after dropping the electron 
masses for the total decay amplitude: 

dr B 



d(dxdydq' dp 2 ± 



G 2 \v ub \ 2 o - 4y 



S{Anfm B y 2 x f 1+ f,± ) 
x\ 4m B x 2 y(l - y) 



-Mo 



2W 2 ' 



Mi + (1 - 2C)M 2 
W 2 



x 1 + (1 - 2C)^ 



-32 M 2 - M 3 x 



using 



Altar,, 



(1 - 2Qm% 
Aui-jr - W 2 



|K^| 2 |p e J 2 



Ami - W 2 



|K±| 2 |Pei| 



V2 8sin 2 (6\y)M^ 



(23) 



Eq. (23) is a simplified form of the decay width in the lim- 
its q' 2 — > and P 2 jml ' B — > 0. In Appendix A we will show 
the full result for finite q' 2 and P 2 . Of the four distribu- 
tion amplitudes g 7rTr , e Tr7r and fj* only the asymp- 
totic form of f Tn is known and there have been attempts 
to model this function and f™ in terms of the instanton 
vacuum [2]. 

A calculation of the other distribution amplitudes is be- 
yond the scope of this article. However, we can see if we 
are able at least to obtain the order of magnitude correct 
if we neglect all contributions except / 7rTr where we posses 
the expression for the asymptotic form and compare the 
semi-leptonic decay B — > iririv with the semi-leptonic de- 
cay B — ► plv. Neglecting e 7 " 1 , K K ' n and g V7r may not be 
so unreasonable, as the first one is twist-3 (i.e. a higher- 
twist contribution) and the other two are connected to 
'polarization' states where we know from the experiences 
of spin physics that the contribution of the quarks is rela- 
tively small. If, in this sense, we only retain the contribu- 
tion from J 71 " 71 ", we obtain a simple expression for for the 
branching ratio B(B — > miiv): 



dB(B -> Tririv) 
d£dxdydq' 2 dp 2 ± 



G 2 \V uh \ 2 



2(4ir) 5 m B r B y 2 x (i + * 
xm%x 2 y(l-y) (M*™*" 2 



(l- q ' 2 ) 
A ) 

i 2 B x 2 y 2 J 



twist -2 = f 
JO 



^e-^0[ C (z, S f)]-^L/r(z) 



where r B is the total B-meson decay width. 



(24) 



4 Numerical estimates 

Now, wc try to make some order of magnitude estimates 
for a comparison between the semi-leptonic decay of a B 
meson into two pions on the one hand and into a p meson 
on the other hand. For / 7r7r one can use the asymptotic 
form given in Ref. [2]: 

ft" as (z, C, W 2 ) = 6z(l - z)(2C - l)F v {W 2 ) . (25) 

Ft T (W 2 ) is the pion form factor in the time- like region, nor- 
malized by Fk(0) = 1. For the pion form factor in the time- 
like region we use the fit from the CMD2-Collaboration 
[16] using the Hidden Local Symmetry (HLS) parameter- 
ization, which is displayed in Fig. 2. The shape of the 
time- like pion form factor is a characteristic superposition 
of the u> and p resonances. 



For the numerics we take Gf = 1.16639 x 10 5 GeV 2 , 
r B x = 1.62 x 1CT 12 s, in B = 5.279 GeV, and Kb = 0.0035 
[17]. The value of V u b is an average value where the error 
assigned to it is of the order of 50%. For the decay width 
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0.8 

W [GeV] 

Fig. 2. Square of the time-like pion form factor (F%) in the 
HLS parameterization. 





M 2 fCJeVl 2 

1V1 [Vjrc V J 
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III}) [Vjrc V J 


left border 


4 


35 


4.7 


central value 


6 


34 


4.8 


right border 


8 


33 


4.9 



Table 1. Sum rule parameters used in the calculation for the 
estimate of the B-decay width. For the calculation of the decay 
width the central values of the sum rule parameters have been 
used, and for the error the maximal deviation to the left and 
right border. 



fs we make use of the corresponding sum rule expression, 
see e.g. Ref. [18]: 



0.0001 



J-2 "*h 

f B = — r ex P 



TO, 



m B - m b 

M 2 



+ 



2M 2 
3 



1 - 



2M 2 



m b {qq) ll 2 =M 2 



(qagGq) „2 =M 2 



8tt 2 



Jral 



sdse m 2 ( 1 



rnt 



(26) 



(qcrgGq)^ = 



-4/A, 



a s (/"o) 

a s (p)\ 2/(3M 



(Mo) J 



(qagGq)^ 



11„ 2 



(w)igcV 2 = -245 MeV 3 
(qagGq) 1 GcV 2 = 0.65 GeV 2 ^! GcV 2 . 



(27) 



The values of the condensates at po = 1 GeV have been 
taken from [14]. In the formula for fs radiative corrections 
are not taken into account because they are not taken into 
account in all the other LCSR calculations presented or 
used here either. For the same reason we use for a s the 
one-loop expression: 



47T 



/?oln( M 2 M 2 ) 



(28) 



using = 0.208 GeV and for the threshold masses m 
the evolution m c = 1.25 GeV and m b = 4.2 GeV [17]. 
The sum rule parameters are chosen as in [14]. They can 
be obtained from Tab. 1 For the calculation of the decay 
width the central values have been used, and for the er- 
ror the maximal deviation to the left and right border. 
For the dependency on the Borel parameter in Fig. 3 we 
consider the differential branching ratio at q' 2 = using 




M 2 [GeV s ] 



Fig. 3. Borel dependency of the differential branching ratio 
B — » -k-kIv. The value for x is fixed at the p-pole in the HLS 
parameterization, i.e. x = 1 — wipHLs/ m s> with m Pi HLS = 
774.57 MeV. 



the central values for the other sum rule parameters and 
integrating p 2 e± e [0, m 2 B x 2 /4], y e [0,1], C & [0,1]. The 
value for x is fixed at the p-polc in the HLS parameteriza- 
tion, i.e. x = 1 — Wp.HLs/ m s> w ith w p ,hls = 774.57 MeV. 
This means x ss 0.97847, so effectively x is close to 1 if 
the invariant mass of the two pions is in the vicinity of 
the p meson pole. It can be seen that the dependency on 
the Borel parameter is rather strong, as we only consider 
the asymptotic form without any higher-twist contribu- 
tions or radiative corrections. Furthermore, for simplicity, 
we have kept the other sum rule parameters fixed. In prin- 
ciple they should vary with the Borel parameter as given 
e.g. in Tab. 1. 

As the p meson decays nearly exclusively into two pions 
there should be a chance to match the branching ratio for 
the semi-leptonic decay B — > nniv with the correspond- 
ing decay B — > piv. More precisely we have to integrate 
the branching ratio over W G [m p — r p , m p + r p ], i.e. the 
p-meson pole and compare: 



dB(B ± -» I/^TT+TT") 



dq 



/2 



dB(B° -» ve-p+) 



q' 2 =0 



4dq 



,2 



(29) 
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Fig. 4. Integration range over the p pole: The shaded area de- 
picts the range of integration in W of the differential branch- 
ing ratio B — » mr£v. 



Here one factor 1/2 accounts for the fact that we have 
to consider a p° wave function instead for a p + and an- 
other factor 1/2 takes into account that we only consider 
charged pions. B(B° — ► p + e~v) at q' 2 = only depends 
on two form factors A\{Q) and A 2 (0) [14]: 



dB(B° -» p+e~v) 



dq' 2 

* q' 2 =0 

A = (m B + m 2 ) 2 — Am B m 
1 



G||Kb| 2 
1927r 3 m|r B 



Vlff 2 



#0 = 



2m 



m B - m p )(m B + m p )A 1 (0) 



A 2 (0) 



(30) 



For the two form factors A\ (0) and A 2 (0) we use the values 
in Ref. [14]: 



Ai(0) = 0.27 ±0.05 
A 2 (0) = 0.28 ±0.05 . 



(31) 



For the integration of the branching ratio B(B — ► -k-kIv) 
over the p pole we take the range W 2 e [(m p — i^) 2 , (m p + 
.Tp) 2 ], see Fig. 4. Hereby we use for the variable transfor- 
mation from x to W 2 the fact that 



W 1 =P 2 =m 



(32) 



We can now compare the the semi-leptonic branching 
ratio B — ► 7T7t^^ with the equivalent £> — > where the 
additional factors are taken into account to make the two 
quantities comparable. Using the central values of the sum 
rule parameters for the calculation and obtaining the error 
from varying the sum rule parameters over the allowed 
range we obtain: 



dB(B° -» p+i-v) 



Adq 



/2 



(I3.8±?fi) x 10- 7 GeV" 



dB{B -» 7T7T^) 



dq 



/2 



(3.0 ±2.7) x 10" 7 GeV" 



g' 2 =0 



(33) 



The consistency in the order of magnitude can be taken 
as a hint that our approach is qualitatively correct. In our 
case the large error results from the fact that we did nei- 
ther take into account all 27rDAs nor the higher-twist cor- 
rections and restricted ourselves to the asymptotic form 
only. The big error in case of the B — > p semi-leptonic 
decay comes from the fact that in the kinematic region we 
consider here, the decay width is the the result of two con- 
tributions that nearly cancel each other, c. f. Eq. (30). The 
number comes from the usual error-analysis, and reveilles 
drastic effects. We should state that this is not a numeri- 
cal estimate of the decay width, but rather a consistency 
check, that we get the order of magnitude correct. The im- 
portant task that remains to be done is a modeling of the 
27rDAs which will then allow for a quantitative prediction 
of the semi-leptonic decay of B mesons into two pions. 



5 Summary and conclusions 

To summarize, we have shown that the semi-leptonic de- 
cay B^ — > Tr + ir~i'l ± can be described in the LCSR for- 
malism using the two-pion distribution amplitudes. Major 
observables, except for / T7r and fj_ n , are the twist-2 dis- 
tribution amplitude j" and the twist-3 amplitude e™. 
When we retain only the twist-2 distribution amplitude 
/ X7r , where the asymptotic form is known, and compare 
the semi-leptonic decay width B — > irirlv with the cor- 
responding decay width B — > piv at q' 2 = 0, we find 
consistency in the order of magnitude. 

I wish to acknowledge A. Khodjamirian and C. Weiss 
for fruitful discussion. 
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x<^ 16 



+M 3 \q' -p, 



M 1 P-p e + M 2 (R-p e ) ll 

(R-q')(R-p e )\\ (P.q')(P.p e 



R 2 

M x P-q' + M 2 R-q' 

(R-q'f (P-q') 2 



pi 



-M 3 I q' 2 - 



R 2 



P 2 



-16 



M 1 P-p e + M 2 (R-p e )\\ 



f (R-q')(R-p e ^ (P-g')(P-p e 
+M 3 I q -p e — — 



A The full formula for the decay width 

In the following we give the formula for the decay width 
B + — > 7r + 7r~ l + vi for finite q' 2 and P 2 . The electron mass 
is neglected in the calculation altogether. The formula can 
be obtained from Eq. (18) by the following substitutions: 



2P-q' 
2p e • q' 
2P-Pe 

2R-q' 

2R- Pe 



2 X-X ,2 

m B x H q 



(m 2 B x-q' 2 )y+- P -^ 



x y 



C -^[m 2 B x 2 + q' 2 ] 

(C-C) 



= (2i? • Pe 



{ 2 /A 

(m B x-q )y - 
+4K ± • p e± 

II + 4K ± • p e± 



x y 



(34) 



In the limit q' — > and i-*0we reproduce the formulas 
Eq. (19). For completeness we add here once more the 
expressions for R 2 , P 2 , and |K^| 2 : 



P 2 = W 2 



m 2 B 



1-x 



R 1 



Ami 



W 



K ± 2 = CCW 2 - ml . 



(35) 



Using the expressions above we can write down the total 
decay width, integrated over the polar angle of p e ± : 

dr B 



d(dxdydq' dpj ± 
G 2 \V ub \ 2 



1 - 



_9 



8(4ir) 5 m B y 2 x ( 1+ pL _\ 



— 32|K _L | |p e j_| M 2 -M 3 



R-q' 



, 2 (R-q') 2 (P-q') 2 



R 2 



■{R-p e )\ + {P-p e f 



P 2 



R 2 



P 2 



p e q 



-2 



, (R-p e )\\(R-q') (P.p e )(P.q') 



R 2 



P 2 



\K ± \ 2 \ Pe± \ 2 (_,2 (P-q') 2 



-4q'- 



R 2 

M\P 2 + MlR 2 + 



q - 



p2 



(„2 ,, 2 R2 \f ,2 (P-q') 2 (R-q') 2 * 



' P 2 



P-q'[(P-p e )(R- q ')-(P-q')(R-p e )\\ 



+P 2 q' 2 



(R.p e ) ll --R.q> 



Mi - M 3 



P-q' 



-[R-q'[(R-p e ) ll (P-q')-(R-q')(P-p e )] 



+R 2 q' 2 



P-Pe-^P-q' 



Mo - M: 



R-g' 



(M 2 W - q' 2 ) 2 



(36) 
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In the form given above the decay width is easy to program 
in a computer code. 



